Metallic Coulomb Blockade Thermometry down to 10 mK and below 
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We present an improved nuclear refrigerator reaching 0.3 mK, aimed at microkelvin nanoelectronic experi- 
ments, and use it to investigate metallic Coulomb blockade thermometers (CBTs) with various resistances R. 
The high-i? devices cool to slightly lower T, consistent with better isolation from the noise environment, and 
exhibit electron-phonon cooling oc and a residual heat-leak of 40 aW. In contrast, the low-i? CBTs display 
cooling with a clearly weaker T-dependence, deviating from the electron-phonon mechanism. The CBTs agree 
excellently with the refrigerator temperature above 20 mK and reach a minimum-T of 7.5 ± 0.2 mK. 



Advancing to ever lower temperatures can open the 
door for the discovery of new physics: for example, 
submillikelvin temperatures in quantum transport ex- 
periments could lead to novel nuclear-spin physics^ in 
nanoscale semiconductor devices'^ or could facilitate the 
study of non-Abelian anyons, Majorana Fermions and 
topological quantum computation in fractional quantum 
Hall samples^i^. However, cooling of nanoscale devices 
below T r-^ 1 mK is a formidable challenge due to poor 
thermal contact as well as microwave and other heating, 
often resulting in device and/or electron temperatures 
raised well above the refrigerator temperature. There- 
fore, significant progress beyond the status quo in both 
cooling techniques and thermometry is necessary. 

One approach to overcome these difficulties uses Ag 
sintersSrJ to thermalize the sample wires2, pioneered by 
the Florida grou p ^"' "K Another approach - pursued by 
our Basel gioup^ - is to use nuclear cooling^"^ on the 
sample wires, with the potential to advance well into the 
microkelvin range. Thermometry in this regime^^^ typ- 
ically faces similar challenges as cooling nanostructures 
and is ideally integrated on-sample. Among numerous 
sensorsi^., Coulomb blockade thermometers^ (CBTs) are 
simple to use and self-calibrating yet offer high accuracy, 
demonstrated down to ^ 20 mK^^. Here, we present an 
improved nuclear refrigerator for cooling nanoelectronic 
samples and use it to investigate CBTs and their mech- 
anisms of cooling. 

We employ a novel scheme for cooling electronic nano- 
structures into the microkelvin regime by thermalizing 
each sample wire directly to its own, separate nuclear 
refrigerator (NR)l^. In this scheme, the sample cools 
efficiently through the highly conducting wires via elec- 
tronic heat conduction, bypassing the phonon degree of 
freedom since it becomes inefficient for cooling at low T. 
A prototype of this refrigerator presented in Ref.[13 has 
been significantly improved in a 2"^^ generation system, 
briefiy outlined below and in Fig.[TJ A network of 21 
parallel NRs is mounted on a rigid tripod intended to 
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FIG. 1. Layout of novel nanosample microkelvin refrigera- 
tor and CBT array. Radiation shields (not drawn) are at- 
tached to the still and cold pate (~ 50 mK). The RC filters 
are 820r2/22nF and 1.2 kJ7 / 4.7 nF. The 21 NR plates are 
0.25 X 3.2 X 9.0 cm'' each, amounting to 64 g Cu per plate. 

minimize vibrational heating. Two separate 9 T magnets 
allow independent control of the NR and sample mag- 
netic field. 

Several stages of thermalization and filtering are pro- 
vided on each sample wire (see Fig.[T]). After 7r-filters 
and thermocoaxiS, each lead passes through a Ag-epoxy 
microwave filter—, followed by an RC filter. Each wire 
then feeds into a Ag-sinter in the mixing chamber (MC), 
emerging as a massive high-conductivity Ag wire. After 
Al heat-switches with fused joints, each lead traverses 
a separate Cu NR via spot-welded contacts, terminat- 
ing in an easily-exchangeable chip-holder plugged into 
Au-plated pins which are spot welded to the Ag wires. 
Therefore, excellent thermal contact (< 50 mil) is pro- 
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FIG. 2. CBT normalized differential conductance g/gr versus 
source-drain dc bias Vsd for various NR temperatures Tcu as 
color-coded, with resulting Tcbt {Sg method, see text) given 
adjacent to each trace. Data from a 67 kQ, a 175 kfi and a 
4.8 MJl CBT is shown. Dashed curves are fits to a model 
(see text). Note lower noise in low-il sensors due to larger 
resulting currents. 

vided between the bonding pads and the parallel network 
of 21 Cu pieces - the microkelvin bath and heart of the 
nuclear refrigerator - while maintaining electrical isola- 
tion of all wires from each other and from ground, as 
required for nanoelectronic measurements. 

The performance of the NRs is evaluated in a series of 
demagnetization runs. The temperature Tcu of the Cu 
pieces is obtained using a standard technique^i^ii^: after 
demagnetization, we apply power on heaters mounted 
on some of the NRs and evaluate the warm-up time- 
dependence Tcu{i) measured with LCMN thermome- 
ters above 2 mK. This allows us to determine both the 
temperature Tcu of the Cu-NRs after demagnetization 
as well as a small field-offset. For each demagnetiza- 
tion run, the NRs are precooled to Ti ^ 12 mK in a 
Bj = 9 T magnetic field and then demagnetized to tem- 
peratures as low as T/ ^ 0.3 mK after the field has been 
slowly ramped down to Bf ^ 0.135 T, giving efficien- 
cies Ti/Tf ^ Bi/Bf > 60%. Reruns showed excellent 
repeatability, allowing us to chart Tcu for various Bj. 
To determine Tcu during the CBT experiments, we use 
the LCMN thermometers above 2 mK, warm-up curves 
at the lowest Bf and in-between, the pre-charted Tcu 
values. 

The network with 21 NRs allows measurements of sev- 
eral CBTs (2-wire each). The CBT devices are Au-wire 
bonded and glued to the Au backplane of the chip carrier 
which is also cooled with a NR. Each CBT consists of 7 
parallel rows of 64 AI/AI2O3 tunnel-junctions in series 
with an area of 2 /im^ fabricated using e-beam lithogra- 
phy and shadow evaporation. The process used allows 
oxidation at elevated temperatures, giving junction re- 
sistances up to IMri^m^. Each island extends into a 
large cooling fin made from Cu, since Cu gives excellent 
electron-phonon coupling. A small B ^ 150 mT is ap- 
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FIG. 3. CBT electron temperature Tcbt versus NR temper- 
ature Tcu for 4.8 MO (open markers) and 67 kfi sensors (filled 
markers, same axes on inset as main figure). Below 10 mK, 
the data is obtained in 3 demagnetization sweeps (blue mark- 
ers) with B = 9T, 5T, 2T, IT and 0.4 T in a typical run, 
ramped at IT/h above IT and 0.5 T/h below. Error bars 
are about the size of the markers. Purple curves are Tcbt 
saturation curves (see text). 

plied perpendicular to the sensor wafer to suppress the 
superconductivity of the Al. The differential conductance 
through a CBT sensor was measured with a standard 
lock-in technique adding a small ac excitation Vac to a 
dc bias Vsd- Note that only 1/64 of the applied voltage 
drops across each junction and the sensor resistance R is 
64/7 times the junction resistance Rj, assuming identical 
junctions. 

We investigated CBTs with various R, see Fig.O Due 
to Coulomb blockade effects, the conductance around 
zero bias Vsd = is suppressed below the large-bias con- 
ductance gx- Both width and depth Sg = 1 — giVsD = 
0)/gT of the conductance dip are related to the CBT 
electron temperature Tcbt- To extract Tcbt, we per- 
form fits (dashed curves) using a numerical model from 
Ref.[l^. We find excellent agreement between model and 
data (see Fig. [2]). Independently, Tcbt can be obtainedi^ 
from the conductance dip Sg — m/6 — u^/60-f m'^/630 with 
u = Ec / {kBTcBT) and charging energy Ec- We first ex- 
tract Ec at high-T assuming Tcu = Tcbt and then use 
this Ec to extract Tcbt from Sg everywhere. While both 
methods produce very similar Tcbt (deviating slightly 
only at the lowest T), the Sg approach makes no a-priori 
assumptions about the cooling mechanism, allowing us 
an unbiased investigation, though now requiring high- 
T calibration against another thermometer (CMN). All 
Tcbt values given here are from the Sg method. 

The thermalization properties of Tcbt of the lowest 
and highest R CBTs are displayed in Fig. [3] for a wide 
range of Tcu from 0.5 mK to lOO mK. As seen, excellent 
agreement is found between Tcbt and Tcu at high tem- 
peratures, as expected. Further, Tcbt is seen to lie well 
above Tcu at the lower temperatures (see Fig. [5] and [3]), 
decoupling fully from Tcu well below 10 mK. We note 
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that Vac was experimentally chosen to avoid self heat- 
ing. Also, the 4.8 Mft sensor reaches lower temperatures 
than the other, lower impedance CBTs, consistent with 
better isolation from the environment, since the power 
dissipated is proportional to V^„^/Rj, with environmen- 
tal noise voltage Venv- 

To model the CBT thermalizationlS, we write down 
the heat flow Qi onto a single island i with electron tem- 
perature Tf. 

(1) 

where Qq is a parasitic heat leak and Vj is the voltage 
drop across the junction, appearing here in the Joule 
heating term. E is the Cu electron-phonon (EP) cou- 
pling constant, = 300 /im'^ the island volume and Tp 
the phonon bath temperature assumed to be equal to 
Tcu- This is well justified by the high thermal conduc- 
tance between the NRs and bonding pads. Note that 
at T <C f K, the sample-to- Au-backplane interface resis- 
tance (Kapitza) is small compared to the EP coupling 
resistances^. Within this model, two cooling mecha- 
nisms are included: Wiedemann-Franz (WF, term) 
and EP cooling. Note the strong dependence of 
the EP term, ultimately rendering WF cooling domi- 
nant at sufficiently low T. Assuming one mechanism and 
simplifying to only one island gives a saturation curve 
TcBT = (T| + T^^J^^P, with a CBT saturation temper- 
ature Ts and an exponent p, corresponding to p = 2 for 
WF-electron cooling and p = 5 for EP cooling. 

We study the mechanism of thermalization by fitting 
the saturation curve first to the 4.8 Mft data. We find 
very good agreement, giving p = 4.9 ± 0.4 (see Fig. [3]), 
indicating that EP coupling presents the dominant cool- 
ing mechanism, limiting Tcbt to 9.2 mK even though 
Tcu — 0.75 mK. Using Qq — ESlT^^j,, a small para- 
sitic heat leak Qq = 40 aW results for each island, with 
E = 2 X 10^ Wm-^K-s from Ref. lS. We speculate that 
Qo could be caused by electrical noise heating such as 
microwave radiation, intrinsic residual heat release from 
materials used or other heat sources. Considering the 
high-i? junctions and correspondingly weak WF cooling, 
it is not surprising that EP coupling is dominant here. 

When analogously examining the low-i? sensors, on the 
other hand, we find p = 3.9 ± 0.4 and Ts = 13.4 mK for 
the 67kn sensor (see inset Fig. [3]), and evenp — 2.7 ±0.2 
and Ts = 6.9 ± 0.1 mK for a 134 kfi sensor (not shown) 
mounted on a conventional dilution refrigerator (base- 
T 5 mK) with improved filtering and chip-holder. Note 
that Ts is the extrapolated T — saturation-T. The 
lowest T measured here was 7.5 ± 0.2 mK. These power- 
laws far below p — 5 indicate that EP cooling is no longer 
dominant but, rather, a more efficient mechanism p < 5 
takes over at the lowest-T in the low-i? sensors. 

In summary, we have demonstrated operation of the 
NRs down to 0.3 mK while the CBTs cool as low as 
7.5 mK. Though the high-i? sensor is obviously cooled 



by EP coupling, the low-i? sensors, interestingly, appear 
to be entering a different cooling regime. However, the 
low-i? sensors have slightly higher Tcbt given the same 
environment, consistent with stronger coupling to the en- 
vironment. The lowest CBT temperatures are limited by 
the parasitic heat leak, which is drained by the cooling 
channels available. 

To further improve the sensor performance, the 
cooling-fin volume can be increased or the heat leak can 
be reduced, potentially using improvements in microwave 
shielding and filtering, e.g. using on-chip capacitors, me- 
tal planes or alternative array designs. Such efforts will 
strongly enhance thermalization if a more efficient cool- 
ing mechanism is indeed present, since otherwise, in the 
EP regime, reducing Qo by 5 orders of magnitude will 
only reduce Tcbt by a factor of ten. 

An alternative avenue based on quantum dot CBTs, 
e.g. in GaAs, might also be rewarding, taking advantage 
of a much larger Ec and level spacing A. The result- 
ing reduced sensitivity to the environment might allow a 
single dot to be used, rather than an array, cooling the 
reservoirs directly via the WF term, rather than through 
a long series of junctions. 
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